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We show that Bogohubov equations in one-dimensional systems with piecewise constant potentials 
can be always solved. In particular, we analyze in detail the case where the condensate wavefunction 
is a real-valued function, and give the explicit expressions for wavefunctions of Bogoliubov excita- 
tions. By means of these solutions, we consider transmission and reflection properties of Bogoliubov 
excitations for two types of potential, namely, a rectangular barrier and a potential step. The results 
yield simple and exact examples of anomalous tunneling effect and quantum evaporation. 
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I. INTRODUCTION 

If the first step toward understanding of a given physi- 
cal system is to elucidate the ground state properties, the 
next step would be the elucidation of the low-energy ex- 
cited states. Indeed, experimentally accessible response 
to an external field is calculated from the information of 
excited states and linear response theory. In many con- 
densed matter systems, low-energy excited states consist 
of not only single particle excitations but also collective 
excitations. As is well known that collective mode in a 
Bose-Einstein condensate(BEC) is well described by Bo- 
goliubov theory^. 

Bose-Einstein condensation in ultracold atomic gases 
has been stimulating many theorists ever since its ex- 
perimental realization^, y, Bogoliubov theory has 
been also well confirmed in these systems by the ob- 
servation of sound propagation and the measurement 
of static structure factor [0, 0|- Another fascinating is- 
sue in ultracold atomic gases is the realization of low- 
dimensional systems. In one-dimension, in particular, a 
variety of physical regimes, from BEC regime to Tonks- 
Girardeau(TG) gas regime, are realized [a^lol. [lOl. [Til, [l^ . 
Theoretically, BEC regime can be attacked by mean- 
field theory and is often referred to as "quasi-one- 
dimensional" [l^ [13] , while TG regime must be treated 
as a strongly correlated system. (Exceptionally, modified 
mean field theory for TG regime has been pro pos ed in 
Ref. [iBl, and its vahdity has been investigated [1^.) In 
this sense, the analysis of Bogoliubov excitations in a one- 
dimensional system, which we will show in the present 
paper, belongs to the studies of quasi-one-dimensional 




FIG. 1: An example of piecewise constant potential. 



'Electronic address: 'takahashi@vortex.c.u-tokyo. ac.jp 



systems. 

In the theoretical works of quasi-one-dimensional sys- 
tems, not only the dynamics of the Bose condensate by 
using Gross-Pitaevskii(GP) equation [13, HI], but also the 
physical properties of Bogoliubov excitations have been 
investigated by using Bogoliubov equations. For exam- 
ple, perfect transmission in the low-energy limit (known 
as "anomalous tunnelin g" ) and related or extended tun- 
nehng problems [l9t l20|, M IM HI, HI, US HI, [13, Hi, 
excitation spectrum and dynamical instability in an op- 
tical lattice [2^, [lOl, and Anderson localization of Bogoli- 
ubov excitations in a random potential[3T|. Recently, 
considering the junction of BECs with different interac- 
tion strengths, an analog of Andreev reflection [s^l and 
Hawking radiation 'ssj are discussed. 

Though the Bogoliubov equations are solved in many 
works numerically or analytically, few exact solutions are 
known. ProbablVjthe solution under a dark soliton or a 
gray soliton 2(J, |33 | is the most important one, in the 
sense that not a few works [22, [2^ [23, [30| have used this 
solution to solve the scattering problem against the po- 
tential barrier with the shape of a delta-function. There 
is no doubt that the delta-function barrier models have 
succeeded in clarifying the fundamental physics of Bogoli- 
ubov excitations. However, several works [231 [25j suggest 
that transmission and reflection properties of low-energy 
Bogoliubov excitations are universal, irrespective of the 
shape of potential. It is therefore highly desirable that 
such universal properties are supported not only by nu- 
merical evidence or zero-width potential (i.e., delta func- 
tion) but also by exactly solved models in the presence 
of barrier with finite width and potential steps. In ad- 
dition, the set of exact solutions will be useful when one 
considers a new physical problem. 

In this paper, we exactly solve Bogoliubov equations 
in a one-dimensional system in the presence of piecewise 
constant potential shown in Fig. [Tj With use of these 
exact solutions, we study the transmission and reflection 
properties of Bogoliubov excitations in two particular ex- 
amples. Stationary GP equation with this kind ofpoten- 
tial shape has been already solved gener ally [13 [11, but 
Bogoliubov equations have not been solved. 

This paper is organized as follows. In Sec. [Til we show 
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the solution of Bogoliubov equations for the real-valued 
condensate function. In Sec. Illli we treat two particu- 
lar examples, i.e., a rectangular barrier and a potential 
step, and solve the transmission and reflection problems. 
In Sec. IIVI we extend the solution when the condensate 
wavefunction is complex- valued. In Sec. |Vl we discuss 
our results and future problems. The conclusion is given 
in Sec. IVll 



II. 



SOLUTIONS FOR THE REAL- VALUED 
CONDENSATE WAVEFUNCTION 

A. Fundamental Equations 



We begin with the following one-dimensional time- 
dependent Gross-Pitaevskii(GP) equation: 



(1) 



In this paper, we consider only the repulsive interaction, 
i.e., g > 0. This equation can be rewritten in a dimen- 
sionless form by introducing the following quantities: 



X 



t 



:t, U{x) 



Uix), 



and ipi^i t) 



(2) 



where ^ is an arbitrary positive constant with dimension 
of length, and often taken to be a healing length. Hence- 
forth, we mainly use these dimensionless quantities, and 
omit bars. Time-dependent GP equation then becomes 



. 9 , 



2dx 



^+U{x)]ij{x,t) + \il;{x,t)\^il;{x,t). 



(3) 



Setting the condensate wavefunction in the form of 

iPix, t) = e-'''*{^'(a;) + [uix)e-''* - v*{x) e'^*] } (4) 

and taking the terms up to first order with respect to 
u{x),v{x), we obtain the stationary GP equation 



L-^ix) = 0, L 



1 d^ 
2dx2 



+ U{x) - fi+\^{x)\^ (5) 



for the condensate wavefunction and Bogoliubov equa- 
tions 



u(x) 
~v{x) 



(6) 



for the wavefunctions of excitations. Needless to say, 
these equations can also be derived by diagonalization 
of mean field Hamiltonian.(e.g., [s^. [37|) 

In this section, we consider only the case where the 



condensate wavefunction is a real- valued function for sim- 
plicity. Discussion on a complex-valued condensate will 
be given in Sec. IIVI Taking account of ^(a:) being real, 
GP equation can be simplified as 



1 d2 



2da;2 

Further, introducing 

S = u + V, G = u — V 

Bogoliubov equations are rewritten as 

HS = eG, 
{H + 2^^)G = eS. 



+ Uix)-fi + ^{xf. (7) 



(8) 

(9) 
(10) 



S and G can be interpreted as phase and density fluctu- 
ations, respectively, as pointed out in Ref. [38i |. 



B. Main Results 

From now, we consider the solution of Bogoliubov 
equations for a piecewise constant potential. It suffices 
for our purpose to consider general solutions for only one 
interval with a constant potential, since the solution valid 
for all regions can be obtained by joining solutions of each 
region smoothly. Thus, we set U{x) = Uq = const. In 
this situation, GP equation ([7]) can be immediately inte- 
grated once: 



(^'')' = *^ + 2(C/o-Ai)*' + CGP, 



(11) 



where the constant of integration Cqp is determined from 
boundary conditions. This equation can be integrated 
once again, and the solution generally becomes Jacobi 
elliptic function. 

By using Cqp in Eq. pip , wc can summarize the so- 
lution of Bogoliubov equations as follows: 

There exists the particular solution such that G is pro- 
portional to S", and S satisfies the following first- order 
differential equation: 



iKdx' 



dS 
K^Jdx 



ie rK'^ 
X V~2~ 



d^- 

da; 



(12) 

5 = 0. 
(13) 



Here the constant K satisfies the following quartic equa- 
tion: 



GP- 



(14) 



Equation determines four possible values of K, 
which is consistent with the number of linearly indepen- 
dent solutions for the original system of equations © 
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and (iTOl) . The solution of Eq. ((T3|) can be always written 
down explicitly even when ^ is an elliptic function, by 
using incomplete elliptic integral of the third kind. See 
Appendix [X] and the example for a rectangular barrier in 
Sec. mi 

We note that this method is applicable even when en- 
ergy e is negative or complex, so Landau instability and 
dynamical instability can be also discussed. 

In the rest of this subsection, we prove the results 
summarized above. Eliminating G or S" from Eqs. © 
and (|10p . we obtain the following fourth-order differen- 
tial equations for S or G: 



HiH + 2^^)G = e^G. 



(15) 
(16) 



We can show that if U{x) — Uq — const, the following 
operator identity holds: 



ax ax 



(17) 



which can be proved by a straightforward calculation 
with use of Eq. ([7]). From this identity and Eq. ([15]), 
we obtain 



ax ax 



(18) 



which means that S' is the solution of ((T6l) . that is, equa- 
tion for G. Therefore, there exists the particular solution 
which satisfies the relation ([T^ . Substituting (fT^ to Eqs. 
^ and pU]) . one obtains 




FIG. 2: Problem of tunneling of Bogoliubov excitations across 
a rectangular barrier. Shaded area represents the condensate 
wavefunction ^{x). 



are solutions of Bogoliubov equations. Here k obeys the 
following equation: 



(23) 



When Ggp = 0, on the other hand, the solutions are 
given by 



' * ^ 2k 



G 



and k satisfies 



i{e + 2(t/o - /i)) ^, 

2ek ^ ' 



(24) 
(25) 

(26) 



HS 



iK 



S', 



{H + 2-^^)8' ^iKeS. 



(19) 
(20) 



Eliminating S'" and S" from HI]), (gOl), and derivative of 
(fT9)) . one obtains Eq. US]). Furthermore, using ((Til), 
(flUl). and derivative of (fT^ . one can obtain Eq. ^4 



C. Special Cases: Solutions Expressed in terms of 
Elementary Functions 

When Ggp = or {Uo — fi)^, ^'(x) becomes an elemen- 
tary function, and the solutions of Bogoliubov equations 
also become elementary functions. In these cases, the 
constant of proportionality between G and S' should be 
taken in a slightly different form rather than (fT2l) to write 
down the solutions neatly. 

When Ggp = {Uq - m)^, = ±(^'2 + Uo- fi) follows 
from Eq. (fTT]) . and 



S = e"=^ ( * ± 



ik 



G = --S' 
2e 



(21) 
(22) 



All possible elementary solutions are summarized in Ta- 
ble |I1 



III. APPLICATIONS 

In this section, we treat two examples, that is, a po- 
tential step U{x) = Uq9{x) and a rectangular barrier 
U{x) = Uo9{a — \x\) {Uo, a > 0), and investigate the 
transmission and reflection properties of Bogoliubov ex- 
citations. 



A. Rectangular Barrier 

As the first example, we consider the tunneling proper- 
ties of Bogoliubov excitations across a rectangular barrier 
U{x) = UoO{a — \x\). See Fig. (21 This problem has been 
first considered by Kagan et aL[21j|. However, Bogoli- 
ubov equations have been solved only numerically. 

Imposing the boundary condition "^{x — > ±oo) = 1, 
one obtains fi = I. (This normalization is always possi- 
ble for any condensate density no at infinity, by taking 
^ = h/y/rrTgriQ in Eq. The solution of GP equation 
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TABLE I: List of solutions expressed in terms of elementary functions. They are classified by the sign of [/q — M and the value 
of Cgp- a, h in the table represent a non-zero real value. If one replaces a; — > a; — i7r/(2a) or a — > io in (ii), one obtains (i) or 
(iv), respectively. Likewise, a; — > a; — 7r/(2a) in (iv) gives (iii), a; — > a; — 7r/(26) in (v) gives (vi), and 6 — > ib in (vi) gives (vii). 
The last (viii) is the limiting case (a, 6 0) of (ii) and (vii). 





Cgp 




*(a;) 


S{x) 


G{x) 


Equation for k 


(i) 




-a^ < 


atanh(aa;) 


ikx j 

^ ( 


a tann(ax) — — 


) 




4 ^ ' 


(ii) 


a coth(aa::) 


ikx J 

^ ( 


acoth(ax) — — 


) 


(iii) 


> 


a tan(ax) 


ikx 

e 


/ / X ifc\ 
' atan(aa:) + — 




4 ^ ' 


(iv) 


acot(aa:) 




zx 


( . X ifc\ 
1 acot(aa:) — —\ 




(v) 







b 

cos (6a;) 


ikx 

e 


fetan(6a:) + i — 


-) 






2 1/7 2 ,2x2 

e = -6 ) 


(vi) 


b 

sin (6a:) 


ikx 

e 


(&cot(&x) 1 2^ 


-) 




(vii) 




b 


ikx J 

' { 


6c 




) 


i(fc' + &')c,' 




sinh(fea;) 


' ^ ^ 2fc 




2ek 


(viii) 








1 

X 




e 


ifcx /I ifc \ 






ifc c' 
2^* 





without node and current is then given by 



of (i) of Table [H They are 



^ J*'"(^) (1^1 <«) 



cn(Q;6a;|m) ' 
^i°"t(a;) = tanh(|.x| - a + tanh"^ 7) 



with 



m = \ — a 



7 



'l + 64(a2-l) 



(27) 

(28) 
(29) 



(30) 



2 + 62(a2 -2)' 
Here the height and width of the barrier are parametrized 



C/o = 1 - 6' 
1 

a — - — cn 

ha 



as 






2 


2 










m 







(31) 
(32) 



Henceforth, we mainly use h and a for parametrization of 
the system rather than a and J7o, because of calculational 
convenience. The case a = 1(<-^ m = 0) is interesting 
since it can be solved elementarily by using (v) of Table 
m However, we do not concentrate on this particular case 
only. 

Next, let us consider the Bogoliubov equations. The 
solution set outside the barrier can be written by means 



S'°"*(x) = e"^'^='(*°"*(x) 



ifci 



■ sgnx 



Gf*(x) 



ifc 



-^[5°"'(^)]'- 
2e ■' 



Here kj {j — 1, 2, 3, 4) are defined as 
ki = V^2(Vl + e2-^, k2 - - 



-ki, 



k. 



iv 2(Vl + + 1), and k^ = -k^ 



(33) 
(34) 

(35) 



The solution inside the barrier is no longer elementary, 
so we must solve Eq. (fT3|) for the elliptic function 4'™(a;). 
Setting the proportionality constant in (fT^ as X = 2e/K, 
one obtains 



5;"(x) - w*i"(x)2 + 



2 , 



4 

462 1 



exp 



-H 



^U&2+K2^a""\^4&2 + K2 



am(a5a::|m) 



G-(x) 
with 



6'/ 2 ON 
-r + -(«'- 2). 



(36) 
(37) 

(38) 
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Here n(n; ip\m) and am(u|m) are incomplete elliptic inte- 
gral of the third kind and Jacobi amplitude, respectively. 
See Appendix [Al nj {j = 1,2,3,4) are roots of the fol- 
lowing equation: 



e'^U^+b'{a'~2)] +b\a'-l). (39) 



More explicitly, we define each Kj 0jS follows: 



4e2 



Kl = 



K3 



62(2-a2) + 2^e2_^64(i_^2)' 

I S2 



62(2 - a2) _ 2^e2 + 54(i_^2) 



K2 = -Ki, 



and K4 = — K3. 



(40) 



Since the sets of the solutions both inside and outside 
the barrier have been prepared, we can now construct the 
solution valid for all regions by joining them. Let us recall 
that our goal is to obtain the transmission amplitude, 
i.e., to construct the solution which has the following 
asymptotic form: 



1 



fc2/(2e) 
fc2/(2e) 



(e"=i^ +re-"=i^) (x ^ -00) 



(x — > +co) 

(41) 



This asymptotic form defines the transmission amplitude 
t and reflection amplitude r. These quantities are well- 
defined, as a result of the constancy of Wronskian: 



W 



, du du* 
dx dx 



. dv dv* 
dx dx 



(42) 



Calculating W from Eq. (gl]), one obtains |tp + |r|2 = 1. 
Physically, the constancy of Wronskian corresponds to 
the conservation of excitation energy [2l| . 

Instead of constructing the solution with the form ((4T|) 
directly, we first construct the even and odd solutions 
free from exponential divergence, because the calcula- 
tion becomes a little easier. Even and odd solutions, 
^^cvcn g-ovcn) {S°'^<^, G°'^'^) , should have the Mlow- 
ing form: 



{x < —a) 



[caFp* -HciF2°"' -^03^4°"* . 
C4(Fl" + F^") + C5(Ff + Ft) {\x\ < a) 

+ 02^2°^' + 03^3°"* [X > a) 



'out 
3 



(43) 



dd 



-^2-^1 - diF°' 
d4(i^r-Fi")-f d5(-F| 

dii^f"* -I- d2i^2°"* + rfa-P?"* 



(x < —a) 
Ft) {\x\<a) 
{x > a) 



(44) 



with F = S and G. From the continuity conditions at 
X = a, the coefficients Ci and di must satisfy 



/r(«)ci+/2™*(«)c2+/3™*(a)c3 

- (ft (a) + ft{a))c, + {ft (a) + ft{a))c,, 
/r(«)rfi+/2"*(«)d2+/3°"*(a)d3 

= iftia) - ft{a))d, + iftia) - ftia))d5, 



(45) 



(46) 



with f = S, S' , G, and G' . The conditions ai x ^ —a 
are automatically satisfied as a result of / being even or 
odd. The solutions of these linear equations can be con- 
veniently written down by the method in Appendix [BI 

Once even and odd non-divergent solutions are 
constructed as the above form, the solution of the 
transmission-reflection problem can be easily derived; 
one has only to make the linear combination of the two 
so that the term proportional to 6"'*^^^ does not exist at 
X +00. It is 



= d 



gcvcn 
2 1 ^cvcn 



C2 



^odd 
^odd 



(47) 



The transmission and reflection amplitudes defined from 
the asymptotic form (j4ip are then given, respectively, by 



t = 



1 2-iki 

2 2 + iki 



1 2-iki 
r = — 



2 2 + iA:i 



d2 
d2 



(48) 
(49) 



Now, we can plot the transmission coefficient |t|2. See 
Fig. [3l It shows perfect transmission in the limit e 0, 
i.e., anomalous tunneling [Tol. [20l. occurs. 

In order to verify the exact solution of Fig. [31 let us 
compare it with reliable established theories in both low- 
energy and high-energy regions. "Lorentzian" in Fig. [3] 
shows the following curve: 



l+72fc2 



da: 



1 

1(2 



1 



(50) 



The author and Kato[39] have shown that this expression 
fits well the transmission coefficient of Bogoliubov exci- 
tations against a high barrier with an arbitrary shape in 
the low-energy region. Indeed, Figure [3l^b) shows that 
this fitting works very well. As the excitation energy e 
becomes high, the transmission coefficient of Bogoliubov 
excitation comes close to that of an ordinary particle 
which obeys Schrodinger equation, as shown in Fig. [Sjc). 
This is because high-energy Bogoliubov quasiparticle has 
almost the same properties as an ordinary particle. 

Figure 2] shows S{x) with very small excitation energy 
e. For small e, G{x) is very small compared to S{x), so it 
is not shown in the figure. It means that for sufflciently 
low energy, u{x) ~ v(x) ~ S{x)/2. The behavior of even 
solution near the barrier is quite similar to that of the 
condensate. Indeed, Bogoliubov equations with e = al- 
ways have the solution {S, G) = (5*, 0) [11], and with the 
use of this solution, perfect transmission has been shown 
for arbitrary shape of potential barrier [2^. 
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FIG. 3: (Color online) (a): Transmission coefficients \t\^ for 
a rectangular barrier with a — 0.71, Uo = 2.96 b — 
0.2, a = 10). Horizontal axis represents not energy e but 
wavenumber k. "Schrodinger" represents the solution for or- 
dinary Schrodinger equation with the same potential barrier. 
"Lorentzian" curve is represented by Eq. (|50p . (b): Close-up 
of low-energy region, (c): Close-up of high-energy region. 



B. Potential Step 

As the second example, let us consider the transmission 
and reflection properties of Bogoliubov excitations for a 
potential step U{x) = Uq9{x). We normalize the conden- 
sate density on the left side as 1, i.e., ^'(x^— oo) = 1, 
which determines the chemical potential as ^ — 1. 

Depending on the height of step [/qj the condensate can 
take the three forms, as shown in Fig.O In the case (a), 
the height of the step is lower than the chemical potential, 
so the condensate on the right side has finite density. In 
the case (c) , the condensate density vanishes in the limit 
X +00 since the step is sufficiently high. The case (b) 
is the intermediate state of (a) and (c), and the conden- 
sate wavefunction shows power-law decay. The constants 
and XL are determined from the continuity of ^'(a;) 
and '^'{x) at x = 0. They are 



_ Jtanh-^( /(2-[/o) /2) (C/q < 1) 
\tanh-i(yT7(2Z7o)) (C/o>l) 



xr 



V2 



Uo 



: COth 



Uo 



[ ^2{Uo ~ 1) 



: sinh 



(51) 

[Uo < 1) 
{Uo = 1) . 
''{2^Uo{Uo-l)) {Uo>l) 

(52) 



2Uo 




5°"(x) 



FIG. 4: (Color online) 5"=™" (a;) and S°'^'^{x) with a = 
2.9, Uo = 0.89 b = 0.4, a = 0.8), and e = 0.04. Upper 
and lower figures show the same plot with different ranges. 
Near the barrier, the even solution is quite similar to the con- 
densate wavefunction. Far from the barrier, they behave as 
simple sine waves. 



(a) Uo<fJ, 



1.5 



-4 -2 

(b) Uo^n 



*(x<0) = tanh(-a;+XL) 
*(x>0) = acoth(a(x+XR)) 



'*(3;<0) = tanh(-a;+a;L) 
1 



*(a;>0) 



-4 -2 

(c) Uo>ll 



X + Xn 



*(a;<0) = tanh(-a;+a;L) 
*(a;>0) ^ 



-4 -2 



sinh(6(a; -f- a^a)) 



FIG. 5: (Color online) Condensate wavefunctions for the po- 
tential step U{x) = Uo9{x). The thick black lines repre- 
sent the potential steps, and the thin blue lines represent 
the condensate wavefunctions. (a) Uo < ji. (b) Uo = ^J.■ 
(c) Uo > M- Note that = 1 in our normalization, and 
a = ^r~Oo , b = ^2(Uo - 1). 



Since the condensate wavefunction is expressed in 
terms of elementary functions, we can use TableUto solve 
the Bogoliubov equations. To avoid the complicated de- 
scription, we prepare several notations before considering 
the transmission and reflection problems. As for the left 
side solutions, referring to (i) of Table [H we define the 
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following: 



Here fcWj = 1,2,3,4) are defined as 



(53) 
(54) 



fcL^iA/2(\/l + e2 + l), andfc:^==-fc^. 



(55) 



The notations for the right side {x > 0) are defined in the 
same way. In the case (a) of Fig. [5l referring to (ii) of 
Table H 

Sf {x) = e"=f " (a coth(a(x + xr)) - ^ j , (56) 



2e 



Gf{x)^--;^[sfix)r, 



kf = V2(v/a4 + e2 _ ^2)^ ^ _fcR 



fc^ = iy^2(Va'* + f"2 + a2), and kf = -kf, 



(57) 



(58) 



with a = \/l — Uq- In the case (b), referring to (viii) of 
Table H 



S 



f(-)=e-f^(— L=-i^), (59) 

(60) 
(61) 



Gf{x) = -^[Sf{x)]', 
kf = iV2i, and fcf = -kf 



In the case (c) , referring to (vii) of Table HI 



5R(a;)^e'"""(bcoth(b(x + XR))-i ^ ' ), (62) 



fcf = ^/2e-b^, kf = -fcf , 



kf^iy/2e + P, and fcf = 



(63) 
(64) 



with b = ^/2{Uq - 1). 

Now, let us construct the solution of transmission and 
reflection problem in the following form: 



(x > 0) 
(65) 




FIG. 6: (Color online) Transmission coefficients for step po- 
tential with step height Uq = 1. Again note that /i = 1 in the 
present normalization. 



The coefficients are determined from the continuity of 
S, S", G, and G' at a; = 0. Since there are four equations 
for five variables, (01,02,04,^1,^3) is determined except 
for an overall factor. It can be explicitly written down 
by using the method in Appendix [b1 The transmission 
coefficient is given by 



T =1- R=l- 



(66) 



Here again, transmission and reflection coefficients are 
well-defined from the constancy of Wronskian(Eq. ([42]) ). 
Note that |di/cip is not a transmission coefficient. 

The transmission coefficients are shown in Fig. IH 
and profiles of wavefunctions of excitations are shown in 
Fig. [7l If the condensate remains finite on the right side, 
Bogoliubov excitations show partial transmission in the 
low energy limit. This situation is shown in Fig.[7I^a) and 
the case Up < 1 of Fig.[6l Such a case is studied in detail in 
Refs. [2^, ll^l in more general systems. The value T(e— > 
0) = 4V1 - Uo/il + VI - Uo)'^ is consistent with these 
earlier works. If the condensate decays on the right side 
because of high step, low-energy excitations show perfect 
refiection, as a free particle obeying Schrodinger equation 
does. It corresponds to Fig. Ulb) and the case Uq > 1 
of Fig. [S] The higher-energy excitations can transmit 
partially, but the component v{x) is suppressed on the 
right side, and only u{x) survives, as shown in Fig.jTjc). 
This means that the quasiparticle is converted into an 
ordinary particle when it is ejected from the condensate 
to the vacuum. Such a phenomenon is called quantum 
evaporation and investigated in detail by several authors 
[13, EH. (However, the present model has no roton-like 
dispersion, unlike the superfluid helium 4.) 



IV. EXTENSION TO THE COMPLEX- VALUED 
CONDENSATE WAVEFUNCTION 

In this section, we discuss the case where the con- 
densate wavefunction takes a complex value. Like the 
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FIG. 7; (Color online) Plot of Re[ti(a;)] and Re[w(x)] of the 
solutions (|65p in various cases, (a) Uo < 1. (b) Uo > 1 and 
e<Uo — l. (c) C/q > 1 and e> C/q — 1. (Note that (M(a;), ^(a;)) is 
multiplied by an overall phase factor so that the amplitudes 
of real and imaginary part become the same.) 



real- valued case in Sec. |TI1 Bogoliubov equations can be 
reduced to a first-order differential equation. Therefore, 
Bogoliubov equations with a piecewise constant potential 
can be always solved irrespective of the form of conden- 
sate wavefunction. 



A. Main Results 

Writing amplitude and phase of the condensate wave- 
function as ^(a;) — A{x) exp[iQ{x)], GP equation be- 
comes 



HA = 0, H 



2da;2 

de 

da; 



2 + U{x) - n + + 



2^4' 



^2' 



(67) 
(68) 



where the second equation is already integrated once, 
and the constant of integration q can be interpreted as 
the condensate current. We re-define S and G as follows: 



G = ue-'^-ve'^. 



(69) 
(70) 



By using these quantities, Bogoliubov equations can be 
rewritten as 



HS- 
{H + 2A^)G- 



i^_d_/G 

Ad^\A 
ig d / 5 
Ad^XA 



eG, 
- eS. 



(71) 
(72) 



Henceforth, we consider the constant potential U{x) — 
Uq = const. We can then integrate GP equation ([67]) 
once, and obtain 

{A'f=A^ + 2A\Uo-ii)-^ + GGP. (73) 

The solution of this equation can be described in terms 
of elliptic functions 35,] . 

With use of Cqp in Eq. ([73)1 . we summarize the solu- 
tion of Bogoliubov equations as follows: 

There exists a particular solution which satisfies the 
relation 



de^ 



(74) 



and S satisfies the following first- order differential equa- 
tion: 



ie [KA^ + Kt - q + {K^ - T)qA-^ ~ g^^-S) 
+ AA' (if 2 + KqA-^ + q{KT - 2q)A-^ - q^rA-^) 
+ SKq + K'^A^ 



S 



q{2q + Kt)A-^ + gVA"* - iqeA-^ A 



S' 
(75) 



with T = K"^ 12 + Uq — ^. The constant K must satisfy 
the following quartic equation: 

= Qif' + C/o - -2Kq~CGP. (76) 

One can easily make sure that if one sets g = in 
Eqs. (HSD-dTSl), they are reduced to Eqs. ([n])-(ITll). Since 
Eq. ((75l) is first-order, it can be solved by the method of 
separation of variables. Thus, it is proved that Bogoli- 
ubov equation with piecewise constant potential is always 
solvable. 

From a practical viewpoint, however, this solution 
seems to be less useful compared to real-valued case, be- 
cause the solution of Eq. ([75|l is difficult to express in 
terms of already known functions. This is in contrast to 
the real- valued case, where the solution can be expressed 
explicitly in terms of incomplete elliptic integral of the 
third kind. 

Before closing this subsection, we sketch the proof of 
the above results briefly. One can show by a direct cal- 
culation that S' + iO'G and G obey the same fourth or- 
der differential equation, even though an operator iden- 
tity corresponding to Eq. (|17p no longer exists in the 
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complex-valued case. Therefore the condition ((74|) fol- 
lows. Eliminating G and high-order derivatives of S from 
Eqs. dm), dZl]), and one obtains Eq. Using 

Eqs. dni), dZSl), (123, and derivative of Eq. dTS]), one ob- 
tains Eq. ((761) . 



B. Check by the solution for gray sohton state 

Though the general solution of Eq. (|75p is difficult to 
express, we can observe its validity by substitution of 
one particular example. Let us consider the gray soliton 
state. (We note that in soliton theory, the term "gray 
soliton" is rarely used, and gray soliton is also called as 
dark soliton.) The condensate wavefunction is given by 




P-^tanh(Wp-^x 



p — ^ ) tanh^ 




P 

Here \q\'^ < holds, q and p represent the supercurrent 
and the density far from soliton, respectively. The solu- 
tion of Bogoliubov equations for gray soliton is given in 
several papers 22, 31, 34]. If it is rewritten in terms of S 
and G instead of u and one obtains 



qk 
2^ 



G 



Akx 



2e 



-A 



where k satisfies 



qk 
2ep 
_ ik 
Ye 

qk 
P 



qk 
iTp 
ik 
1 
qk 
P 



fc2 

y 



qk 



A' 



p-A^ 



1 



From expressions ((79)) and ((80)) . one can show 

G. 



5' - i ( - - ^ - ^ 
* p A^ 



(79) 



(80) 



(81) 



(82) 



Furthermore, one can make sure by direct substitution 
that Eq. ([TS]) with K = (2e/A:) - [q/ p) has the solution 
of the form d79l). 



V. DISCUSSION AND FUTURE WORKS 

A. Relation to soliton theory and extension to 
spinor BEC 

When we solve the Bogoliubov equations. Equa- 
tions (|12p and ((7^ are the most important ansatz. We 



have derived these conditions by showing that G and S" 
(in the complex-valued case, G and 5" -I- i6'G ) obey 
the same differential equation. A more sophisticated un- 
derstanding of these ansatz comes from Lax pair of soli- 
ton theory. In Ref. ^3^], it is shown that the solution of 
Bogoliubov equations can be represented by the square 
of the solution of linear problem of Lax pair, called as 
"squared Jost function" . Using this fact and one of lin- 
ear equations (Eqs. (7) and (9) of their paper), one can 
derive the equation which is equivalent to our ansatz, 
namely, Eqs. dH]) and ((7^ . 

The method of squared Jost function suggests the pos- 
sibility of the extension to spinor BEC. In Refs. (i^. |43|. 
113], it is shown that spin-1 Gross-Pitaevskii equation be- 
comes integrable when the strengths of interaction sat- 
isfy a certain condition. The integrability of this sys- 
tem originates from that of matrix nonlinear Schrodinger 
equation [4^. We conjecture that Bogoliubov equa- 
tions of this integrable spin-1 BEC system can be also 
solved by means of squared Jost function, as well as 
scalar BEC. Investigation of this conjecture and appli- 
cation to the tunneling problem of collective modes are 
future works. 



B. Instability of soliton train states under the 
perturbative potential 

In Sec. IIII[ we have mainly applied the exact solutions 
to the transmission and reflection problems of collective 
excitations. Another important application of our solu- 
tions is a stability study of soliton train states [13, Is^ . 
Stationary CP equation in a uniform and infinite sys- 
tem has the solution of the form ^'(x) cx sn(aa;]m). Un- 
der this snoidal condensate wavefunction, the solutions of 
Bogoliubov equations can be obtained by using Eq. ((13)) 
and Appendix [X| From the explicit expressions of ex- 
act solutions, one can derive the following facts: (a)For 
positive energy < e < ec, where Cc is a certain posi- 
tive value determined from chemical potential, two of the 
four linearly-independent solutions have "positive norm" , 
i.e., lu(a;)]^ — lf(a;)l^ > 0, and the other two have "neg- 
ative norm", i.e., lM(a:)l'^ — lw(a;)l^ < 0. (Even though 
the solution in an infinite system cannot be normal- 
ized, the sign of lM(a:)l^ — |'y(2^)|^ is meaningful.) (b)For 
e > ec, two of the four solutions have "positive norm", 
and the other two are unphysical solutions which diverge 
exponentially. Since Bogoliubov equations always have 
the solution pair (e,u,v) and (— e,w*,M*), the above (a) 
indicates that the soliton train state possesses Landau 
instability |47|. It should be noted that ec goes to in the 
limit of sn{ax\m) — > tanh(aa;), i.e., one dark soliton limit. 
Therefore, one dark soliton state has no Landau instabil- 
ity. (This is clear from the explicit expression in (i) of 
Table HI) Thus, the existence of negative energy and pos- 
itive norm eigenstate is a property specific to soliton train 
states. We further note that positive norm state and neg- 
ative norm state degenerates. As shown in Ref. 48], this 
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degeneracy is a necessary (not sufficient) condition for 
emergence of complex eigenvalue. Therefore, when per- 
turbative potential becomes sufficiently strong, one can 
expect that the dynamical instability would appear, as 
demonstrated in a double well system[49j. This rough 
discussion on infinite system would be modified in a fi- 
nite system, because eigenstates are discretized. In order 
to obtain more precise results, a more quantitative study 
for a finite system with perturbative potential is needed. 



C. Supercurrent state through a potential step 



tunneling effect and quantum evaporation. 
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In the studies on physical origin of anomalous 
tunneling (lol. [20l . [2l| . it was important to reduce the 
physical property of low-energy excitations to that of 
the condensate wavefunction. Indeed, the perfect trans- 
mission has been proved for generic potential by using 
the coincidence between the condensate wavefunction 
and the wavefunctions of excitations in the low-energy 
limit ^2^. Furthermore, the similarity between the super- 
fiuidity of the condensate and the perfect transmission 
of low-energy excitations has been proposed [2^. These 
works suggest that the perfect transmission of excitations 
can be regarded as an extended concept of the superflu- 
idity of the condensate. 

In the presence of potential step, as shown in Refs. [2^, 
[27| and Sec. lIIII of the present paper, the excitations show 
partial transmission at zero-energy. (See Fig. [S] again.) 
Recently, the present author has found that there exists 
a supercurrent state through a potential step without re- 
flection, whose density profile is quite similar to Fig.[5fa). 
We conjecture that Bogoliubov excitations show partial 
transmission in the low-energy limit in this supercur- 
rent state, as the case without supercurrent does. If so, 
it gives a counterexample to the above earlier work(26j 
which relates the superfluidity of the condensate and the 
perfect transmission of the excitation. That is to say, the 
existence of the supercurrent state of the condensate does 
not necessarily mean the perfect transmission of the ex- 
citation. Thus, this problem may give a chance to recon- 
sider the physical interpretation of anomalous tunneling 
effect. We will report on this issue elsewhere. 



VI. CONCLUSION 

In this paper, we have shown that one-dimensional 
Bogoliubov equations with piecewise constant potentials 
can be reduced to a first-order linear differential equation, 
and therefore, it can be always solved by the method of 
separation of variables. Particularly, when the conden- 
sate wavefunction is a real-valued function, the solution 
of Bogoliubov equations can be expressed explicitly in 
terms of the elliptic integral. Using these solutions, we 
have solved transmission and refiection problems of ex- 
citations for a rectangular barrier and a potential step. 
Our results provide new exact examples of anomalous 



APPENDIX A: FORMULA FOR INCOMPLETE 
ELLIPTIC INTEGRAL OF THE THIRD KIND 

We only show the formula necessary in this paper. 
More information on elliptic integrals and elliptic func- 
tions is available, e.g., in Refs. [50l[5l|. 

Jacobi amplitude is defined as the inverse function of 
incomplete elliptic integral of the first kind: 



am ^{(p\m) 



d0 



'o \/l -m sin^ 9 
Jacobi elliptic functions are then defined as 

sn(u|m) — sin(am(u|m)), 
cn(u|m) — cos(am(u|m)), 
d 

dn(u|m) = — — am(M|m). 
ou 



(Al) 



(A2) 
(A3) 

(A4) 



Incomplete elliptic integral of the third kind is defined as 



n{n;ip\m) = — 
Jo (1 



{1 — n sin'^ 6)\/l — msin^ 9 
Setting if = am(u|TO) and 9 — am(z|m), one obtains 

dz 



(A5) 



n(n; am(u|m)|m) 



1 — nsn2(z|m) 



(A6) 



Since any squares of Jacobi elliptic functions are related 
to sn^, one can always solve Eq. by means of this 

formula. 



APPENDIX B: SOLUTION OF n - 1 
HOMOGENEOUS LINEAR EQUATIONS WITH n 
UNKNOWNS 

Consider four homogeneous linear equations with five 
unknowns: 



^aii ai2 ai3 ai4 aisN 

021 0^22 ^23 024 ^25 

031 O32 033 O34 035 

\a4i 042 043 044 045^ 



X3 
\X5J 



(Bl) 
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Note that the coefficient matrix is not square. Obvi- 
ously, except for an overall factor, {xi,X2,X3,X4,X5) is 
determined uniquely. The solution is given by 



Thus, 



Aj is the solution of Eq. (jBl 



A, 



where A,- is defined as 



A, = 



aij+2 

fl2,i+l «2,i+2 0,2, j+3 0,2,j+i 

+ l a3,i+2 a3j+3 ^3^+4 

O-i.j + l a4j'+2 a4,j+3 a4,j+4 



(B2) 



(B3) 



Here the matrix indices are considered by mod 5, i.e., 
j + n is replaced by j + n — 5 if greater than 5. It can 
be easily proved by cofactor expansion of the following 
trivially zero determinant: 



Oil 0.i2 tti^ QiQ 

ail ai2 ai3 an ois 

121 022 023 ^24 025 

031 032 033 034 035 

a4i a42 043 044 045 



^a,,Aj=0 (i = 1,2,3,4). 

(B4) 



When one of Xj, say X5, is fixed to a certain value, this 
formula is reduced to well-known Cramer's rule. Gener- 
alization of the above result to n — 1 homogeneous linear 
equations with n unknowns is straightforward. 

This formula is convenient to write down the solution 
of linear equations in a compact form. For example, ci 
and C2 in Eq. (|45p are given, respectively, by 



Cl 



C2 



^out / 
^out 
^out / 

COUt 
QOUt / 

^3 
^out 

^out / 



^3 

COUt / 

•-'3 

^out 


-Gi" 




cin 

-ST' 
-GT 


Q\n 

0^ 

oin / 

-Gi" 


^out / 


-Gi"' 


- G'f 


-cr 


— LT4 


-S'T 


-G'" 


-Sf 

cin / 

-Gf 


-s'T 
-s'T' 

— 


^out / 
^out 


-Gi"' 




-Gf' 


— LZ4 


^out / 



(B5) 



(B6) 



where the arguments of all matrix elements are a, which 
are omitted to save space, di and c?2 in Eq. can be 
obtained in the same way. By using these expressions, the 
transmission coefficient (j48p can be written down com- 
pactly and explicitly. 
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